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INTRODUCTION 

Critical gas turbine engine hot section components such as blades, vanes, and 
combustor liners tend to develop minute cracks during the early stages of 
operation. These cracks may then grow under conditions of fatigue and creep to 
critical size. Current methods of predicting growth rates or critical crack sizes 
are Inadequate, which leaves only two extreme courses of action. The first is to 
take an optimistic view with the attendant risk of an excessive number of service 
failures. The second is to take a pessimistic view and accept an excessive number 
of "rejections for cause" at considerable expense in parts and downtime. Clearly 
Is is very desirable to develop reliable methods of predicting crack growth rates 
and critical crack sizes. 

To develop such methods, it is necessary to relate the processes that control 
crack growth in the immediate vicinity of the crack tip to parameters that can be 
calculated from remote quantities, such as forces, stresses, or displacements. The 
most likely parameters appear to be certain path-independent (PI) Integrals, 
several of which have already been proposed for application to high temperature 
Inelastic problems. A thorough analytical and experimental evaluation of these 
parameters needs to be made which would include elevated temperature isothermal and 
thermo-mechanical fatigue, both with and without thermal gradients. 

Investigations of fatigue crack growth under elastic-plastic condition should 
consider the impact of crack closure on the appropriate crack growth model. 
Analytically, this requires the use of gap elements in a nonlinear finite element 
code to predict closure loads. Such predictions must be verified experimentally 
through detailed measurements; the best method for measuring crack closure has not 
been established in previous studies. 

It is the purpose of this contract (NAS3-23940) to determine the ability of 
currently available Pi-integrals to correlate fatigue crack propagation under 
conditions that simulate the engine combustor liner environment. The utility of 
advanced fracture mechanics measurements will also be evaluated and determined 
during the course of the program. These goals are to be accomplished through a 
nine task, combined experimental and analytical program. To date, an appropriate 
specimen design, a crack displacement measurement method, and boundary 
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condition simulation in the computational model of the specimen has been achieved. 
Computational verification of the path-independence of the proposed Integrals has 
been demonstrated for Isothermal and thermal gradient cases. Also, the 
experimental testing and data acquisition is continuing. Tensile and cyclic tests 
were run at several strain-rates so that an appropriate constitutive model could be 
developed. The experimental data include cyclic crack growth tests under 
isothermal, thermo-mechanical, and thermal gradient conditions. 


EXPERIMENTAL PROGRAM 

Alloy 718, a y-y" nickel-base superalloy, has been selected as the analog 
material for this program because over the temperature range from 800 to 1200°F, it 
shows very large changes in creep behavior. This permits the use of Alloy 718 to 
simulate the behavior of combustor liner materials while still performing 
experiments at a relatively low temperature. Tensile, creep and cyclic 
constitutive tests have been performed over the temperature range from 70 to 
1200°F. Even though this material can experience large amounts of creep 
deformation at the upper end of this temperature range, the tensile and cyclic 
tests showed little evidence of strain rate sensitivity on constitutive response. 
Figure 1 shows the resultant tensile stress strain curves determined for 70, 800, 
and 1200°F. Tensile curves were also determined at 900, 1000 and 1100°F. These 
curves were used during the analysis of the temperature gradient tests which are 
described later In this paper. 

The crack growth tests have been performed using a single edge notch (SEN) 
specimen with buttonhead grips. The width and thickness of the gage section is 0.4 
and 0.1 inch, respectively. Detailed descriptions of the specimen, testing 
procedure and finite element analyses have been described previously [1]. The SEN 
tests were performed in a strain control mode with the experimental setup shown 
schematically in Figure 2. The controlling extensometer was mounted at the center 
of the 0.4 Inch wide surface of the specimen. The other two displacement gages, 
one to monitor crack mouth opening displacement and one to monitor the back face 
deflection, were also used. The controlling and back face extensometers had a gage 
length of 0.5 inch and the crack mouth gage had a gage length of 0.03 inch. The 
crack length was monitored using a DC potential drop technique. The hysteresis 
loops from these tests [1] were open, indicative of large net section non-linear 
deformation (i.e., plasticity and creep). Figure 3 shows the data from four 
1000°F, A c =® crack growth rate tests - single tests at strain ranges of 0.5 
and 1.7% and duplicate tests at 1.15% strain range. The initial EDM notch depth 
was approximately 0.01 inch. Figure 3a shows the strong influence of strain range 
on the crack growth behavior. The results from the duplicate 1.15% tests are very 
similar. During these strain control tests, the maximum loads diminished as the 
crack grew. Figures 3b and 3c show the 1000°F crack growth rates (da/dN) plotted 
as a function of AK and K max , respectively. These data show that these tests 
cover a range of conditions which can not be described by K, which is proportional 
to the crack tip stress field in linear elastic fracture mechanics. Comparison of 
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Figures 3b and 3c also shows that the crack growth rates are better described by 
K max than aK - This stJ 99 ests that when the final analysis of these tests is 
completed, the adjustments for crack closure will describe the behavior over the 
entire range of crack length. 

Thermal gradient tests have been performed to evaluate the Pi-integrals under 
non-isothermal conditions. The thermal gradient was established using a 
combination of induction heating and forced air cooling. At the front or cracked 
surface, the crack mouth opening extensometer prevented the use of cooling air so 
the crack was grown from the high temperature to the low temperature. Cooling air 
was applied to the back face of the specimen above and below the extensometer with 
cooling tubes and to the gage section with air forced through the cooling passages 
in the extensometer probes. This gradient technique was originally developed on a 
specimen which was monitored with fifteen thermocouples. Three sets of five 
thermocouples were attached to the specimen along the crack plane and at planes 
0.25 inch above and below the crack plane. The two latter locations correspond to 
the positions where the arms of the two large gage length extensometers contact the 
specimen. On all three planes, a thermocouple was located at each edge of the 
specimen and at four equally spaced Intervals. The presence of the thermocouples 
prevented the attachment of the three extensometers, so the temperatures were 
monitored along the plane of the crack in another specimen using an optical 
pyrometer attached to a traveling microscope. The results of the temperature 
measurements are shown in the Figure 4. The line connecting the "X" symbols 
represent the optical pyrometry measurements. These results show a small amount of 
temperature difference from specimen-to-speclmen and along the gage length. This 
variation is within an acceptable range. 


COMPUTATIONAL PROGRAM 

Major computational accomplishment was the implementation, verification, and 
applications of the selected Pi-integrals to uniform and non-uniform thermal 
gradients, and thermo-mechanical loadings for fracture mechanics analyses involving 
nonlinear material stress/strain behavior. A critical review [2] of the available 
J x -integrals conducted under this program revealed that Pi-integrals proposed by 
Ainsworth (Jq), Blackburn (J*), Kishimoto (J), and Atluri (Tp* and Tp) are 
suitable for nonlinear thermo-mechanical response. The relative advantages and 
limitations of the various J x -integrals for thermal gradient problems are 
discussed in-depth in Reference [2]. The traditional Rice J-integral becomes 
path-dependent and loses its physical significance for thermo-mechanical loadings 
due to presence of the mechanical (rather than total) strain energy in its 
formulation and, also, due to the absence of thermal strains contribution and 
non-homogeneous (temperature dependent) material properties. 

Verification of the implemented thermo-mechanical J x -integrals was carried 
out by prescribing a linear temperature gradient through the width of the SEN 
(single edge notch) crack specimen which was gradually loaded with far-field 
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uniform normal stress. The crack length to specimen width ratio (a/W) Is 0.25 and 
the specimen length to width ratio (L/W) Is 1 . 25 . The specimen edge at crack face 
was kept at 1200°F temperature and the back-face edge had a temperature of 
1000°F. All properties of the material were kept temperature dependent and, as 
such, they varied through the specimen width. Therefore, the material properties 
are non-homogeneous along the direction of thermal gradient. 

Nonlinear, thermo-mechanical, finite element (FEM) analyses of the SEN specimen 
were carried out for several values of the applied uniform stress, resulting In 
thermo-elastic to thermo-elastic/plastic response. Figure 5 shows the far-fleld 
normal displacement (Uy) variation which Is linear In nature through the specimen 
width for all the load cases considered. A significant amount of plastic yielding 
is taking place ahead of the crack for the larger values of the applied stress as 
seen In Figure 6 which shows the effective stress contours in the specimen. The 
FEM analyses results were then used In the post-processor Pi-Integral computer code 
to determine the various J x -1ntegrals for several integration paths around the 
crack tip. Table 1 shows a typical variation of the various J x -1ntegrals along 
four paths spanning extreme plasticity to pure elastic regions near the crack. ,It 
could be seen that the Rice J-integral values are path dependent for the 
thermo-mechanical load case. The other five J x -1ntegrals shown in Table 1 are 
path-independent. A few of these thermo-mechanical Pi-integrals are plotted in 
Figure 7 for various applied normal stress values. These integrals show excellent 
agreement with isothermal Elastic -Plastic Fracture Mechanics (EPFM) Handbook [3] 
and Tada/Paris [4] analytical solutions for lower values of applied stresses. For 
higher values of applied stress, the thermo-mechanical Pi-Integrals attain larger 
values as compared to analytical isothermal J-integral. These results demonstrate 
implementation and verification of the thermo-mechanical Pi-integrals considered 
here. 

To minimize the finite boundary effects, the specimen L/W ratio used in 
verification analysis was 1.25. For the actual specimen, however, the L/W ratio is 
only 0.78. Therefore, in order to investigate the influence of finite boundary on 
J x -integrals, the same linear thermal gradient analysis was carried for the 
actual specimen. Typical results as shown In Figure 8 indicate a maximum deviation 
of 3% on Pi-integrals for the prescribed applied stresses. 

The actual thermal gradient developed in the SEN specimen gage-section can be 
approximated as trilinear relationship, as shown in Figure 9. It is approximately 
constant at 1200°F for the first 0.175 inch along the specimen width; it then 
linearly d^ops to 1050°F for the next 0.14 inch of the width; and, finally, it 
varies linearly to 900°F in the remaining 0.08 inch of the width. Stress 
malysis of the specimen with the measured temperature dependent material 
properties and the prescribed trilinear thermal gradient was performed. An 
interesting result was observed for the pure thermal gradient load case (zero 
applied load/deflection). For an a/W ratio of 0.25, a crack mouth opening 
;i '.placement of 0.00015 inch was predicted. The corresponding far-field normal 
displacement for pure thermal expansion load case is shown in Figure 10 for two 
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cross sections - one at y = 0.31 Inch (end of gage length), and the other at y = 
0.25 Inch (location of extensometer probes). These free-thermal -expansion normal 
displacements vary nonllnearly, through the specimen width, due to a combination of 
the temperature dependent material properties and the trlllnear thermal gradient. 

It has been shown previously that the far-fleld mechanical displacement In the 
specimen varies linearly along the specimen width at the end of gage-section. 

Figure 11 shows the normal stress variation ahead of the crack-tip for the pure 
thermal gradient load case. The normal stress has a value of 20 ksl near the 
crack-tip and it drops sharply to -8 ksl at x = 0.175 inch and then gradually 
increases to 12 ksl by the end of the specimen width. This normal thermal stress 
is self-equilibrating in nature since there is no mechanical load applied. This 
fact is further elaborated In Figure 12 which shows the contour plots of the normal 
stress In the entire gage-section of the specimen. It could be seen that at the 
top surface (end of the gage-length) the normal stress is zero along entire width 
of the specimen. The whole specimen Is In elastic state of stress with a maximum 
effective stress of 35 ksl. Since there exists a crack-tip stress field In the 
specimen for the measuredlthermal gradient load case, the various J x -integrals 
were determined and found to be path-independent. The average value of these 
PI -integrals was 1.16 lbs/inch, which Is equivalent to 5 ksi/in thermal-Kj 
value. 

Stress analyses were further carried out by prescribing uniform normal stresses 
superimposed on the actual specimen with the measured thermal gradient. Figure 13 
shows the variation of the thermo-mechanical J x -1ntegrals with applied stress. 

It could be seen that the values of these integrals are very close to each other 
for lower values of the applied stress with very little plasticity in the 
specimen. When the level of plasticity increases, then the various Pi-integrals 
have somewhat different values but remain path-independent. 


CONCLUDING REMARKS 

The analyses conducted for thermo-mechanical Pl-integrals have been very 
successful and compare very well with the available results. Work is now in 
progress for computational simulation of the measured load-displacement plots for 
the thermo -mechani ca 1 response of the specimen. In addition, computations are also 
to be performed for reversed plastic loading to simulate crack-closure. These 
results will be compared with crack growth experimental measurements to identify 
path-independent integrals which can predict cyclic crack growth behavior under 
isothermal, TMF , and thermo-mechanical loading conditions. 
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Figure 1: Tensile Curves for Alloy 718 at 70, 

800 and 1200°F. 



Figure 2: Schematic Drawing of SEN Specimen 

Test Method 
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FIGURE g: ’h«rmo-Mtchan!c*1 Klshtnoto tntaaral In SEN So*c<n«n 
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